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We continue the investigation of thermodynamical properties of the BPS Skyrme model. 

In particular, we analytically compute the baryon chemical potential both in the full field 
theory and in a mean-field approximation. In the full field theory case, we find that the 
baryon chemical potential is always exactly proportional to the baryon density, for arbitrary 
solutions. 

We further find that, in the mean-field approximation, the BPS Skyrme model approaches 
the Walecka model in the limit of high density - their thermodynamical functions as well as 
the equation of state agree in this limit. This fact allows to read off some properties of the 
w-meson from the BPS Skyrme action, even though the latter model is entirely based on the 
(pionic) SU(2) Skyrme field. On the other hand, at low densities, at the order of the usual 
nuclear matter density, the equations of state of the two models are no longer universal, 
such that a comparison depends on some model details. Still, also the BPS Skyrme model 
gives rise to nuclear saturation in this regime, leading, in fact, to an exact balance between 
repulsive and attractive forces. 

The perfect fluid aspects of the BPS Skyrme model, which, together with its BPS proper¬ 
ties, form the base of our results, are shown to be in close formal analogy with the Eulerian 
formulation of relativistic fluid dynamics. Within this analogy, the BPS Skyrme model, in 
general, corresponds to a non-barotropic perfect fluid. 

PACS numbers: 


I. INTRODUCTION 

The derivation of properties of baryonic matter, especially beyond the nuclear matter density, 
is still one of the most challenging problems in current strong interaction physics. Since neither 
perturbative nor lattice computations apply in this regime, one is forced to use an effective model 
approach where both the field content and the form of the action are postulated from general 
considerations (symmetries, low energy degrees of freedom) rather than derived from the underlying 
fundamental quantum theory, i.e., QCD. 

One of the most popular and successful effective theories is the Skyrme model framework [1], 
where the fields are reduced to low energy effective chiral fields (pions and, optionally, higher 
mesons). Baryons, on the other hand, are not introduced as independent degrees of freedom but 
appear, instead, as collective excitations in this mesonic matter, that is, as topological solitons 
with an identification between the baryon charge and the topological degree. 

Very promising results of the original Skyrme proposal when applied to the baryon sector El-S 
(and also to some light nuclei) have to be contrasted with its problems in the modeling of higher 
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nuclei and (infinite) nuclear matter. There are two main reasons for this fact. Firstly, the nuclear 
binding energies found in the Skyrme model are too large. Secondly, skyrmions with high baryon 
charge behave as crystals, that is, form a lattice built out of smaller charge substructures [9], which 
is in contrast to the liquid-type behavior of nuclear matter. Recently, three possible ways to cure 
the binding energy problem have been proposed: one may include a dominating sextic term in the 
model ( the near-BPS Skyrme model |lf)j . [] !’. [12]). add (infinitely) many vector mesons ( the BPS 
Skyrme vector meson model ps], m ), or include a ’’repulsive” potential (the lightly bound model 

ra¬ 
in particular, the BPS Skyrme model (the BPS restriction of the near-BPS Skyrme model) al¬ 
ready provides quite accurate binding energies of the most abundant higher nuclei (after taking into 
account the semiclassical rotational and iso-rotational corrections as well as the Coulomb interac¬ 
tion and a small isospin breaking EE]). It also leads to a perfect fluid energy-momentum tensor 
with SDiff symmetries. Therefore, this model plays the role of a field theoretical realisation of the 
liquid droplet model of atomic nuclei. Its near-BPS generalization, thus, seems to be a natural 
candidate for an effective model of nuclear matter within the Skyrme framework. Furthermore, the 
fact that the BPS Skyrme model has the energy-momentum tensor of a perfect fluid (without any 
averaging or mean-field approximation) led recently to a noticeable progress in the understanding 
of thermodynamic properties of skyrmions and nuclear matter as a skyrmionic medium |l7j. With 
our model, it was possible to find an exact equation of state relating the local energy density and 
the pressure which, as the energy density is generically spatially dependent (local spatial fluctua¬ 
tions), allows to study skyrmionic matter beyond the mean-field approximation. This was crucial 
for a better understanding of the long-standing issue of the too high compression modulus (too low 
compressibility) in the Skyrme model [TTJ. The reason for these achievements is a rather unique 
property of the BPS Skyrme model, namely the equivalence of the microscopic (in terms of observ¬ 
ables related to the original effective chiral fields) and macroscopic (in terms of thermodynamical 
functions and variables) thermodynamical descriptions. 

The sextic term, which is just the square of baryonic current, has been included in many effective 
models [18j (see recent m, and may be induced by the topological WZW coupling with the vector 
meson [20]-[25j. It is known to significantly improve quantitative predictions of the model in 
the baryon sector. 

The fact that the BPS Skyrme model provides the rare possibility to study nuclear matter 
beyond the mean-held limit has been recently employed in a neutron star context, where the BPS 
Skyrme model was coupled to gravity |26j (the original Skyrme model was used for the description 
of neutron stars in [27]). It has been found that such skyrmionic stars have masses and radii in very 
good agreement with current data for neutron stars, but, what is perhaps even more important, 
their properties do change if one performs the full held theory and gravity computation instead of 
the usual mean-held approximation known as the TOV approach. 

Beyond this progress, there is, however, another thermodynamical function (besides the energy 
(density), pressure and volume) which is extremely important but poorly understood in Skyrme 
type models. Namely, the baryon chemical potential, which is crucial if one wants to consider 
the model as a realistic description of cold and dense nuclear matter. It governs the stability of 
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phases, possible phase transitions as well as the coexistence of different phases. Therefore, for 
any realistic application of gravitating BPS Skyrmions to neutron stars (which means taking into 
account /3-equilibrium, existence of skin and crust or possible quark core - hybrid stars) one has to 
know this thermodynamical function. 

It is the main aim of the present work to fill this gap in our knowledge of the thermodynamical 
properties of the BPS Skyrme model. 

We shall find that it is again the geometrical nature of the BPS Skyrme model which allows 
to find the baryon chemical potential both in a mean-field approach and in an exact calculation, 
providing further evidence for the importance of geometric models of nuclear matter. Furthermore, 
using these results we are able to compute the mass of a BPS skyrmion in a skyrmionic medium, 
which is of high phenomenological importance as it may be related to the in-medium masses of 
baryons. 


II. PROPERTIES OF THE BPS SKYRME MODEL 

The BPS Skyrme model is defined by the following Lagrange density (we use the metric con¬ 
vention (+, -, -)) 

C BPS = jC 6 + jC 0 = —(24tt 2 ) 2 A 6 B^ - X 0 U. (II.l) 

It consists of the sextic term i.e., the baryon current squared 

B" = ^e^Tr L v L p L a , L p = U^8 P U (II.2) 

and a non-derivative part, that is, a potential U(JJ). In the full near BPS Skyrme model, this 
BPS action is supplemented by the usual Skyrme (perturbative) model. However, in order to keep 
the binding energies on an acceptable level, the BPS part must give the leading contributions to 
masses, while the non-BPS part should enter as a rather small addition. Hence, for the description 
of static properties already the BPS Skyrme model should give a rather accurate approximation. 
For convenience, we redefine the coupling constants as A6 = A 2 /(24) 2 and Ao = v 2 ■ 

Starting from now, we restrict considerations to the static case. The BPS property of the model 
means that one can reduce the static field equations to a first order (Bogomolny) equation 

Avr 2 H 0 = ±vVti (II.3) 

whose solutions saturate the following topological bound 

Ebps > 2tt' 2 Xu\B\(VU), (Vti) = dnVU (II.4) 

27T z J g3 

where (VU) is the average value of the square root of the potential on the target space. Further, 
B = J d 3 xB° is the baryon number (topological degree) of the Skyrme field. In fact, this equation 
can be analytically solved and possesses infinitely many SDiff related solutions in each topological 
sector. 
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As already mentioned, the energy-momentum tensor has a perfect fluid form and, for static 
configurations, reads 

T 00 = e, T ij = - P5 ij , (II.5) 

where the energy density e and pressure P are 

e = A 2 t t a BI + v 2 U, P = A 2 t t a BI - v 2 U. (II.6) 

Hence, the Bogomolny equation is equivalent to the zero-pressure condition, [28] . Furthermore, 
the equation with a non-zero pressure gives rise to (non-equilibrium, one parameter) solutions of 
the full equations of motion, where from the conservation law for the energy-momentum tensor 
it follows that P must be constant. Furthermore, the energy density - pressure equation of state 
(EoS) is not an algebraic equation but depends on the field component (off-shell) 

e = P + 2 v 2 U (II. 7) 

or, after inserting an exact solution, on the space coordinates (on-shell) 

e = e (P 5 £) (II.8) 


where the particular dependence of the energy density on x depends both on the form of the 
potential and on the particular solution. This is a very peculiar property of the BPS Skyrme 
model: one can find an EoS (for the skyrmionic matter in the thermodynamical equilibrium) 
without any mean-field approximation i.e., generically, with a non-flat energy density distribution. 
In this sense, the BPS Skyrme model is a unique field theoretical tool which goes far beyond usual 
EFT. 

The fact that, generically, one deals with spatially non-constant energy densities poses, at first 
glance, a problem if one wants to compare the results of the BPS model with typical EoS obtained 
for other EFT (in a mean-field approximation). However, it is possible to derive a mean-field 
equation of state in the BPS Skyrme model, as well. 

Obviously, one may define an algebraic EoS which connects the total energy E and pressure: 
E = E(P) or geometrical volume V and pressure: V = V(P). Interestingly, by means of the 
Bogomolny equation, one can rewrite the total static energy and the geometric soliton volume as 
target space integrals independently of any particular solution and finds 

E(P) = 2ir\v\B\E, V(P) = 2 tt\B\^V (II.9) 

where 


or 


E= fries in 2 g 2 ^±4 , V{P) 
Jo ffU + P 



d£ sin 2 £ 


1 

Vu + P 



2U + P \ 
s/U + P/ ’ 



( 11 . 10 ) 


Vu + P 


(ILll) 
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and P = P/v 2 . Here, in (II. 10) we used the usual Skyrme field parametrization 


U = cos£ + i sin£ n ■ t , n = (sin© cos 4>, sin 0 sin <h, cos 0) 


( 11 . 12 ) 


where t are Pauli matrices and n is a unit vector. Further we assumed that the potential depends 
on the Skyrme field only via the profile function £. It is sometimes useful to consider the axially 
symmetric ansatz for a skyrmion with baryon number B, where 

f = £(r), 0 = 0, <f> = Bp (11.13) 

and (r, 9, ip) are spherical polar coordinates. Then the Bogomolny equation takes the form 


\B\X 
2 r 2 


sin 


ftr = -WU + P. 


(11.14) 


We want to emphasize, however, that all thermodynamic functions and variables as well as the 
thermodynamic relations between them are completely independent of any ansatz. They follow 
directly from the BPS equations (or, more generally, once integrated static field equations) and 
are the same for all solutions. They represent, therefore, properties of the model itself, and not of 
particular solutions. 

The bulk observables E(P) and V(P) are given by the average values (on target space) of certain 
functions of the potential ’’shifted” by the pressure. They obey the thermodynamic relation 


P = - d l 

’ dV 


(11.15) 


and, therefore, coincide with the corresponding thermodynamical functions, which means that 
for the BPS Skyrme model the microscopic (field theoretical) description is equivalent to the 
macroscopic approach (by thermodynamical functions). Despite the fact that, generically, we do 
not have an algebraic density-pressure EoS (which is an important fact as it allows to go beyond 
the mean-field approximation), one can perform such a limit and define an average energy density 


E 

~ £= v 


(11.16) 


which obviously possesses an algebraic relation to the pressure (or volume). The fact that in the 
model we can compare mean-field with non-mean-held computations has been used recently for 
neutron stars [29). In particular, we studied how their properties are modified by going beyond 
mean-held and taking into account the spatial dependence of the energy density of nuclear matter. 
In general, these two formulations (MF and exact) allow to investigate which predictions stem from 
the model itself, and which are related just to the mean held approximation. 

Another important quantity is the particle number density which here is just the baryon charge 
density 


Pb = Bo- 


(11.17) 


Again, usually it has a non-constant (spatially dependent) form. An obvious proposition for an 
average particle (baryon) density is 

PB = y, B = f d 3 xB 0 , 


( 11 . 18 ) 
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as the total number of particles (the baryon charge) in a given volume V is B. 

Although the energy density and pressure (or the particle density and pressure) are not related 
by an algebraic equations of state, there is an algebraic relation which connects all three of these 
local quantities, namely 

e + P = 2A 2 7T 4 p|. (11.19) 

This equation will play a prominent role in the present work. 


III. BARYON CHEMICAL POTENTIAL 


A. Definition and properties 


The standard thermodynamical definition of the chemical potential is provided by the following 
relation with other thermodynamical functions 


From (11.19) it follows that 


e + P = pp B 


p = 2A 2 7 T 4 Pb- 


(III.l) 


(III.2) 


Hence, the baryon chemical potential is proportional to the baryon charge density and, similarly 
to this quantity, it is a spatially non-trivial function. Let us remark that this is an off-shell, i.e., 
solution independent result. 


Again, as for the energy density and particle number density, we may define a mean-field 
chemical potential. In general, it reads p = ( dF/dN ) where F is the free energy, and N is the 
particle number. In our zero temperature case we have F = E and, further, the particle number 
N is the baryon number B. We, therefore, get 


M = 



(III.3) 


Of course, p must necessarily obey the same thermodynamical relation (for the averaged quantities) 


£ + P = PPB- 


(HL4) 


This definition should be understood as follows. We consider a skyrmion in equilibrium, i.e., a 
solitonic solution with a given baryon charge Bo which occupies a volume Vo and has an energy 
Eo . Such a solution is a solution of the BPS equation and, therefore, corresponds to zero pressure 
P = 0. Now, we want to check how the energy of a skyrmion varies if we increase its topological 
charge to B = Bo + n but keep the volume fixed. Of course, this cannot be done in a smooth way, 
since the baryon charge is a topological, conserved quantity. Nonetheless, we may find a solution 
of the pressure equation which has increased topological charge B and occupies the same volume 
V 0 . The price we pay for that is the appearance of a non-zero pressure, which depends on the 
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’additional’ baryon charge, P = P(n). Then we get two equations (here Vo does not change, i.e., 
V 0 (n,P) = V o (n = 0,P = 0)) 


E(n) = 2 -k\v(Bq + n) / sin" £ 


A 7 71 

Vq = 2ir(Bo + n) - / sin 2 £ 


2 ^2 U + P 


Vu + P 

1 


y/u + p' 


Using the definition (III.3) of p and the fact that the volume remains constant we find 

2 U + P 


TV 


p = 2ir\v— 
r 2 


\fu 


+ P 


1 


Vu + P, 


(III.5) 

(HI-6) 

(III-7) 


This can be simplified to the following formula which expresses the chemical potential as a function 
of the pressure 


7r 


p = AtiXu— 
2 


u + 4 


(HI.8) 


On the other hand, equation (III.7) leads to the proper thermodynamical relation 

P = flp B ~ £ 


(III.9) 

which proves the consistency of our definition of the baryon chemical potential with standard 
thermodynamics. It is instructive to compare this formula with (11.19). If we integrate (11.19) over 
the soliton domain we get 


PV + E = 2vr 4 A z J <Pxp%. 

Hence, the MF baryon chemical potential can be found in the following, alternative form 


1 


p = -^27t 4 A 2 / d 6 xp z B = 27t 4 A" 


_ o 4 \ 2 


/ d 3 xpl 


(III.10) 


(III.ll) 


/ (PxpB 

which, after applying the Bogomolny equation, agrees with the target space average derived above. 
Obviously, p / 2\ 2 it 4: pb, except for a very special case, that is, the step-function potential, for 
which all local quantities coincide with the corresponding average (MF) ones. 

Furthermore, as 

\ -i 


B 1 v2 
PB = TT = 


1 


V 0 2vr A 7T \ ^Ju + P/v 2 


we get 


£ = p 


1 v 2 


1 


-l 


- P 


(III.12) 


(III.13) 


2vrAvr \- s /U + P/v 2 t 

which allows to express the energy density as a function of the chemical potential. Comparing now 


formulae (III.8) and (III. 13) we re-derive the well-known relation 

dp ' 


dP 


1 

V PB 


(III. 14) 
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B. High pressure limit 

In general, the baryon chemical potential - both the exact one and its mean-field counterpart - is 
a complicated function of the particle number (baryon density), whose detailed form is governed by 
the used potential. The same is true for other thermodynamical quantities as well as the equation 
of state, which, too, depend on the potential. However, one can observe that at high pressure 
(equivalently energy or particle density) the model reveals a universal behaviour. Namely, it tends 
to the BPS Skyrme theory with the step-function potential. Indeed, for P > i/ 2 we get at leading 
order 


E = tt 2 A By/P (III.15) 

V = tt 2 B-^=. (III. 16) 

Hence, 

£ = P + Boq, pB = —(III.17) 

7r-A 

and 

p = 2ir A \ 2 p B (III.18) 

where is a bag constant at infinite pressure, see [29 ]. The fact that, asymptotically, the baryon 
chemical potential grows linearly with the baryon charge density (particle density) is generic for 
the BPS model and is not affected by a particular form of the potential. 


C. Low pressure limit 

At vanishing pressure the MF chemical potential is always equal to the equilibrium energy Eq 
divided by the topological charge 

Ao = ^ (III.19) 

while the exact chemical potential is, up to a multiplicative constant, the baryon density at equi¬ 
librium. This happens even in the non-compacton case i.e., when solitons are infinitely extended 
and the geometrical volume is infinite, which leads to zero average energy and particle density. 

For compact skyrmions, the model realizes a liquid-gas phase transition. Indeed, at zero pres¬ 
sure and in the given (equilibrium) volume Vo, we may have solutions with a smaller amount of 
topological charge by just removing some compactons from the given volume Vo. All these solu¬ 
tions are stable and form a collection of compact solitons, surrounded by empty space. As they 
are BPS solutions (P = 0), their total energy is exactly proportional to the topological charge and, 
therefore, p = pQ. Obviously, the corresponding MF energy density and MF charge density will 
tend to 0 as the topological charge decreases. 
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D. Thermodynamic relations for specific potentials 

1. The step function potential: U = 0 (Tr (1 — U)) 

We start with a very special case - the step-function potential. This is a unique choice for 
the potential in the BPS Skyrme model which results in a constant energy density and particle 
number density (baryon charge density). Therefore, all local quantities completely agree with their 
mean-field (averaged) counterparts. The step-function potential reads 

U = 0 (Tr (1 - U)). (III.20) 


The energy and pressure are constant 


E{P) = l^±J=^, V{P) = l ' 


2 y/l + P' 


2 y/l +P 


(III.21) 


and 


E(P) = 2tt\u\B\E(P), V{P) = 2it-\B\V(P). 


(III.22) 


Let us compute the baryon chemical potential from its definition (III.3). As before, we assume that 


we start in the equilibrium (where P = 0) and then we add more particles (increase the baryon 
charge) keeping the volume constant. Let 


E 0 = E(P = 0) = 27T 2 Ai'|L>o|, Po = V(P = 0) = tt 2 ^\B 0 \ 


Then, after changing from Bq to B = Bq + n we get 


(III.23) 


p 0 = tt 2 X B- 1 


1 + ™ 


(III.24) 


2 2 + ^#) 
E = tt 2 \uB- 


1 + ^1 

jyZ 


(III.25) 


where P is due to the higher, than at the equilibrium (B = Bo), topological number. Both E and 
P are functions of B. Simplifying this expression we find 


* (S).-(s;,. 


However, from the constant volume condition we get 


( dP 


2tt 4 A 2 


\dBjy V/ 


B 


and 


(III.26) 


(III.27) 


[i = 27t 4 A 2 ^- = 2ir 4 \ 2 pB- 
bo 


(III.28) 
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v 2 *e / v 2 *P 



FIG. 1: Average energy density (continuous line) and pressure (dashed line) as functions of the chemical 
potential for the step function potential. 


Similarly one can find 


D 4 \ 2 -2 2 1 2 2 

P = n\p B v = ^ 4x2 p v 

(III.29) 

4 \ 2 -2 , 2 1 2,2 

e = ir A Pb + v = +u 

(III.30) 

(see Fig. 1). Obviously, the relation 


P + E= pp B 

(III.31) 

holds. Moreover we get the usual EoS 


£ = P + 2u 2 . 

(III.32) 


As we see, the baryon chemical potential is always proportional to the particle density, not only 
asymptotically as has been proven in generality. Furthermore, the local and global chemical po¬ 
tential equations are exactly the same and p = p. 

2. No potential: U = 0 

It is also possible to find the chemical potential for the BPS Skyrrne model without any potential. 
Of course, at the equilibrium (no pressure) there are no stable soliton solutions, but for any non-zero 
P skyrmions do exist. Then, 

E(B) = tt 2 XBVp, V = ir 2 \B-^= (III.33) 

where we assume that we start with a given non-zero pressure solution with a fixed topological 
charge Bq and then increase B keeping the volume constant. Hence, 

E(B ) = 7r 4 A 2 -^- =>- p = 27t 4 A 2 — 


(III.34) 
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H = 2it 4 \ 2 pb 


(III.35) 


Moreover, as 


P = £ = -K 4 X 2 p 2 B 


(III.36) 


we get that 


P £ 47t 4 A 2 ^ 


(III.37) 


This case is quite similar to the step-function potential as the energy density and baryon density 
are again spatially constant. All quantities can be obtained from the former case by a simple v —> 0 
limit. 


3. Cubic potential: U = g(£ — \ sin 2^) 

Let us now consider a more nontrivial situation, that is, a potential which leads to non-constant 
energy density. A simple example can be provided by the cubic (in the sense of the approach to 
the vacuum) potential U = \ sin2£). This potential belongs to the so-called BPS potentials 

and provides exact and particularly simple solutions of the non-zero pressure integrals. Then, 


V = 2 


or explicitly 


Hence 


+ P - y/p] , E = - 


- \/P 


E=- -P \ p + 


2 ttJ- + P-PV 


+ 7 T\ P + 


£ = y r 


v f + ^ 

/ 7 T , P 
2 ' V 


/ P V* 
7? 


(III.38) 


(III.39) 


(III.40) 


Then, performing similar computations as before, we find how the energy and pressure vary if the 
baryon number is changed 


P = v 


.pB 2 / B 2 \ 2 0 tt pi 


1 _ Po) = ,,2 77 Pb f i Pop 
B* 8h 2 ol 


(III.41) 


/ B 2 1 Bi B 2 Y\ 

+ bT 4 Bp 1 - Bp j 

Here po,s is the average baryon density at equilibrium. Therefore, 


Be B 3 


Po,b p B 


(III.42) 


(III.43) 
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FIG. 2: Left figure: Average density as a function of the MF chemical potential for the cubic potential 
U = — \ sin 2£). Here v 2 = A = 1. Right figure: zoom-in close to saturation density. 


P 



FIG. 3: Pressure as a function of the MF chemical potential for the cubic potential U = |(£ — 4sin2£). 
Here u 2 = A = 1 


and 


£ = 




lBl 

4 B 2 





1 Ab (, . P 2 b 
4 Pb V P°o,b 


(see Figs. 2,3). Asymptotically, for sufficiently high chemical potential we get 


and 


or 


M = 


7T 5 / 2 

V2 


Xu 


B 

B~o 






u 


2 


P = 


1 -2 
4^ 4 A 2/i ’ 


8 = 


1 -2 
47t 4 A 2 ^ ’ 


(III.44) 


(III.45) 


(III.46) 


(III.47) 
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Hence, for high value of the chemical potential we as always obtain the linear EoS 

P = e (III.48) 

where the subleading constant B^ has been omitted. 


4 - Non-compacton potential: U = |(£ — ^ sin 2£) 2 


Another potential we want to discuss is a potential with a sextic approach to the vacuum. It 
means that BPS skyrmions are no longer compactons but usual, infinitely extended solitons. It 
results in an infinite geometrical volume at the equilibrium. Then, the average energy density goes 
to zero at the equilibrium. Hence, similar as in the non-potential case, we are forced to close the 
skyrmionic medium in a given, finite volume by the application of a nonzero external pressure. 
Then, we find 


V = 27 t—B arsinh——, 
v 2 y[P 


E = ir \vB\l -1- — 

v 4 


P 2 


7T 


Thus, in terms of the baryon number 

P = 


vrV 1 


7T 


4 sinh 2 ^F 2 


, E = —A vB cotlr 


v V 
2 -kX B 


or using the average particle density 


P = 


ir 2 v 2 


4 sinh' 


2ir\pB 

Therefore, the MF chemical potential reads 

_3 


1 7 P 

2 —-—, E = —XvBcoth. 


2itXpb 


- n \ 
a = —Azz 

2 


coth 


+ 


2-kXpb J ' 2-nXps sinh 2 j 


For high particle density we find the usual linear relation 


p = 2 it a X 2 pb, Pb -t oo. 


(III.49) 


(III.50) 


(III.51) 


(III.52) 


(III.53) 


Here it is also possible to derive the MF chemical potential for vanishing particle density 

^ i/ l 

p = — Xv + 7 T 2 zz 2 —e pb : p B —> 0. (III.54) 

2 p B 

Since the average mean-field energy density goes to zero at the equilibrium (P = 0), this case, at 
least from the mean-field perspective, does not look like a bag type model. The saturation density 
is simply zero. 
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5. The Skyrme potential: = \ Tr (1 — U) 

Finally, we want to present the result for the usual (pion mass) Skyrme potential 

Z4=^ = ^Tr (!-[/) = 2sin 2 ! 


(III.55) 


Then, the MF baryon chemical potential as a function of the average baryon density is implicitly 
given by the following formulae 


^ = AnXis 


15i/ 4 


/ P f A o o. 

\ 2z/2 1 


\ 2zy2 1 

J 2 + -2 ( (A^ + 2v 2 P + P 2 )E 

2V 2 + P 

- P(v 2 + P)K 

2u 2 + P 


(III.56) 


Pb = 


3 v 


8vrA 


2 +r[(^ + p)E 


2v 2 

2 u 2 +P 


- PK 


2v 2 

2 u 2 +P 


(III.57) 


Here, E and K are the complete elliptic integrals of the first and second kind, respectively. 


IV. IN-MEDIUM SKYRMIONS 


Using the framework presented in this paper, it is also possible to obtain energies (masses) of 
skyrmions in the skyrmionic medium. This is of high importance, as it allows to find in-medium 
masses of baryons (nucleons). In order to find how the energy of a skyrmion with a given topological 
charge Bq varies if it is immersed in the skyrmionic medium, we consider the following integral 


2<h 2 U + P 


Eb 0 (p) = 2tt\uBq / d£sin 2 g-_ 

Jo VU + P 


(IV.l) 


This is the energy of charge Bq skyrmion under the external pressure P. This pressure is induced by 
the additional baryon charge n injected into the equilibrium solution (volume Vo, pressure P = 0) 
without changing the volume 

A 


V 0 = 2ir(B 0 + n) — I dg sin 2 g 

u Jo VU + P 


(IV.2) 


If compared with equation (III.5), the energy differs only by the change of the overall multiplicative 


factor Bo + n into Bq- This corresponds to the fact that now we are interested in the energy of the 
charge Bq baryon surrounded by a skyrmionic medium with the additional charge n, which leads 
to a non-zero pressure. Of course, this pressure is the same in the medium and in the original 


(now compressed) skyrmion. Therefore, the energy (IV.l) is, in fact, the in-medium energy of the 


charge Bq skyrmion. Using these formulas and the expression for the baryon chemical potential, it 
is possible to express the energy of a skyrmion as a function of the mean-field chemical potential. 
As an example, we consider the step function potential. Then, one can show that 


Eb 0 {b.) = n 2 \v(B 0 + n) ( 1 + 


R2 

n o 


(B 0 + n) 2 


(IV.3) 
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which, after taking into account the formula for the baryon chemical potential, can be rewritten as 


p , n B o (, 4tt 4 AV 




(IV.4) 


(remember p = p for the step function potential). This is valid for p > po i.e., above equilibrium. 
For p = no, i.e., from the vacuum value until the equilibrium, the in-medium mass is always the 
same and equal to the equilibrium mass. To show this, let us assume that we start with a collection 
of charge one skyrmions with the total baryon charge Bq, which in equilibrium occupy the volume 
Vo- Now, we take away charge one skyrmions one after the other. Obviously, due to the contact 
form of the interaction and the BPS nature of the solutions, each removed unit skyrmion has the 
same energy (mass). However, the density of the medium, which is now a gas of BPS skyrmions 
in the fixed volume Vo, decreases. 

Repeating the computation from the previous section, we can conclude that for any potential 
asymptotically the in-medium energy of a skyrmion (baryon) always behaves as 


e b 0 (v) = at p ->• oo. 


(IV.5) 


In general (arbitrary potential), the in-medium energy of charge one baryon reads 


E b = i = M- =-) (IV.6) 

1 + n p B 

where the second part vanishes at the equilibrium (P = 0) and tends to ^p at asymptotically 
large densities. Let us underline that this result is beyond the mean-field (constant energy density) 
approximation. 

In the full near-BPS Skyrrne model or in the BPS Skyrme model with semiclassical contributions 
included, it is reasonable to expect a modification of the obtained in-medium mass dependence. 
Indeed, due to non zero binding energies one may expect a skyrmion mass which increases with 
decreasing medium density at small densities. At sufficiently high density, the universal relation 
e b 0 ~ should again be valid. 

In an analogous manner, we may compute the in-medium size of a BPS skyrmion. In general, 
the volume Vo occupied by the original charge Bq soliton at equilibrium is reduced by adding the 
additional topological charge n i.e., by increasing the medium density 

Vbo = Vo = — -> 27T _^ Bo at p ->• oo, (IV.7) 

Bo + n pb p 

and the radius of the compacton reads 


Bbq 



3tt 2 A 2 

2 



at p -» oo. 


(IV.8) 


Observe that, due to the thermodynamical properties of the BPS Skyrme model, the global quan¬ 
tities (the mean-field energy density, the mean-held particle density, the mean-held chemical po¬ 
tential) of a skyrmion and the surrounding medium always coincide. On the other hand, their local 
counterparts differ. 
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V. PHYSICAL IMPLICATIONS 

So far, we have presented some thermodynamical properties of the BPS Skyrme model from a 
more theoretical point of view. In this section, we want to study their effect on physical properties 
of baryons and nuclei if described by BPS skyrmions. Furthermore, we shall use our results to 
relate the BPS Skyrme model to effective theories including c o and a mesons and show how their 
properties are hidden in our skyrmionic (pionic) action, with an exact balance between attractive 
and repulsive channels. 


A. The Walecka model and BPS skyrmions 


The Walecka model of nuclear matter consists of nucleons (neutron and proton spinors) which 
interact with the scalar a meson and the vector u meson m 

Cw = C-N + + Cint (V.l) 

where 

C N = $ - mil 11 7°) ip (V. 2 ) 
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(V.3) 


C-int = gaipcn/j + (v.4) 

Here u Fl/ = d F oo u — A self-interaction term for the scalar meson may also be included. The 

mean-field approximation means that one computes the partition function (in the thermodynamical 
limit) 

Z = J V^Vi/iVaVu ef Cw (V.5) 

in the limit where the bosonic fields are approximated by their constant condensate values a, ui o- 
Then all derivative dependent terms disappear and the interactions are simplified to a mesonic 
background field seen by the nucleons. This means that we arrive at a free fermion model with 
shifted parameters, 

£w = 4> “ m *N + ^*7°) ( V - 6 ) 

where 


m* N = m N - g a a , //*=//- 0 . 


(V.7) 


One should remember that the baryon chemical potential (which enters in all thermodynamical 
relations) is still fi. However, the effective chemical potential fi* sets the Fermi energy of the 
’’effective” free fermions 


E* f = n* = 



(V.8) 
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Here we assume the zero temperature limit. The baryon density reads m 


2 k 3 F 

PB= 3^ 

and is related to the co meson vacuum value by 


9uj 

wo = — pPb- 

mi 


(V.9) 


(V.10) 


An interesting observation is that in the limit of large density, i.e., kp —> oo, the equation of state 
in the Walecka model exactly coincides with the MF EoS derived for the BPS Skyrme model at 
pressures P S> B^ (which in practice means P v 2 ). Namely 


P = £. 


(V.ll) 


Moreover, in the Walecka model this limit reads as 


e _ 1 91 2 

6 “ o 2 PE' 
2 mi 


(V.12) 


Comparing this with the universal relation between average energy density and baryon density in 
the MF BPS Skyrme model we find that 

4 A 2 1,^, ( V .l3) 

2 ml K ’ 

Another important observation can be made if we analyse the large density limit for the effective 
chemical potential. Then we find 


kp —> oo 


p* = kp and p* 


1/3 

Pb 


(V-14) 


Hence, at kp —> oo 

p — p + fl'wWo — p H- kPb ~ — ?~Pb- (V-15) 

mi mi 

This formula has an exact counterpart in the BPS Skyrme model in the MF approach and in the 
high density limit 


p = 27t 4 A 2 pb- 


(V.16) 


Comparing the last two expression, we find independently again the relation (V.13). It is a striking 
fact that in the exact (non-mean field) microscopic thermodynamics in the BPS Skyrme model this 
formula is valid at any pressure (density). Indeed, as we know 


p = 27t 4 A 2 pb- 


(V.17) 


The conclusion is that the BPS Skyrme model and the Walecka model are equivalent in the high 
density regime, as far as thermodynamical properties are concerned and a mean-field approximation 
is made. The obtained relation between the BPS Skyrme and Walecka model parameters allows 




18 


us to read off the ui meson coupling constant in the BPS Skyrme model. For ra u = 738 MeV and 
g 2 /47r = 10 — 12, which is the empirical value of the cu-nucleon coupling, we get the acceptable 
values for the Skyrme model parameter A 2 = 9 — 11 MeV fm 3 . It may be also compared with an 
upper bound for the coupling constant = 25.4 [32]. Then, A 2 < 47 MeV fm 3 . Let us remark 
that all previously used potentials led to A 2 slightly bigger than the optimal value, but significantly 
below the upper bound [Tbl [26l [29] . Undoubtedly, this bound on the BPS Skyrme model parameter 
may help to constrain the potential part of the action. 

The reason why these two models possess the same thermodynamical properties at high density 
limit originates from the fact that in this regime the Walecka model is dominated by the vector 
meson sector. From the interaction point of view, the vector meson couples to the baryon current 
and effectively the dominating part of the Lagrangian is the baryon current squared. But this 
is exactly the derivative part of the BPS Skyrme model, although the latter is expressed in a 
topological manner and not by nucleon spinors. We comment that the same effect shows up in a 
quark bag model with a vector interaction [331 . In any case, it is a rather interesting observation 
that it is possible to read-off the ratio between the u meson mass and its coupling constant to the 
nucleon in the BPS Skyrme model, even though there are no obvious u meson degrees of freedom 
in its action. Instead, we only have the usual pionic fields. The u meson is hidden in the form of 
the BPS action. We remark that the sextic term of the BPS Skyrme model may also be obtained 
from a Skyrme type model with the oj meson included explicitly [32], as the leading contribution 
in a derivative expansion [20]. This leads to the plausible conjecture that soliton solutions in a 
Skyrme model with uj mesons should be quite similar to soliton solutions in a Skyrme model with 
the sextic term (the baryon current squared) included. That this is indeed the case has been 
demonstrated recently numerically in the baby Skyrme model in one lower dimension [34] . where 
the baryon current squared is identical to the Skyrme term. 

At lower densities, the EoS of both models get much more complicated, and the BPS Skyrme 
model EoS now depends on the potential, such that a direct comparison is no longer obvious. In 
addition, at low densities the field theoretical picture of the BPS Skyrme model is probably more 
adequate than the mean-field approximation. For example, we know that the constant energy 
density approximation is not the proper one if one wants to derive the compressibility of skyrmionic 
matter at equilibrium un¬ 
it is interesting to note that the local chemical potential is exactly linear in the baryon density 
for the whole range of the pressure, not just in the asymptotic regime, as happens in the Walecka 
model. Besides, one may perhaps expect that there is a relation between the equation of state in 
the BPS model and a particular form of potential for the a meson field. 

B. In-medium mass of nucleons 

Now we want to identify nucleons with charge one skyrmions, in order to see what results we 
get from our in-medium skyrmion mass calculations for the corresponding in-medium masses of 
nucleons. As we have shown, the mass of a charge one skyrmion is density independent as long as 
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the medium density does not exceed the saturation density, i.e., the energy density at equilibrium. 
Above this density, it starts to grow with an asymptotically linear dependence M ~ This 
behaviour agrees with the in-medium mass behaviour recently found using a holographic approach 
[35], where the mass of nucleons decreases from its vacuum value until one reaches the saturation 
point, where it begins to grow. If we want to compare with further approaches, the following 
subtlety must be taken into account. By construction, our in-medium nucleon masses are always 
the total in-medium skyrmion (rest) energies per baryon number. This is not the case in several 
approaches (e.g., in the Walecka model, see Eq. (V.7), or in [36]), where the in-medium mass 
of the nucleon is exclusively induced by in-medium changes of certain coupling constants, which 
may even lead to an in-medium reduction of the nucleon mass above nuclear saturation. In these 
approaches, further contributions to the total energy stem from in-medium modified nucleon- 
nucleon interactions, which may give significant contributions to the total in-medium energy per 
baryon number. Obviously, our results must always be compared with these total energies per 
baryon number. 

Within our framework, the reduction of the nucleon mass between the vacuum and saturation 
density expected on general physical grounds, may probably be obtained by an extension of the 
model to its near-BPS generalization by taking into account the perturbative part of the Skyrme 
theory (pionic perturbative part) and/or by the semiclassical corrections. 

Let us also remark that the in-medium properties of skyrmions in the present work, in contrast to 
the in-medium modified Skyrme Lagrangian m, have been obtained within the BPS Skyrme model 
without introducing any medium-dependent constants. Such a medium modified BPS Skyrme 
model would have the following form 

£bps = —7t 4 A 2 {fj,,x)BcrB a — v 2 (n,x)U. (V.18) 

Using our results one may try, however, to express the coupling constants A and P as functions 
of the medium density /r (and, probably, of the coordinates x ) and fit to the correct in-medium 
mass dependence. Then one may check how good this approximation is by comparing with other 
thermodynamical properties. At this stage of research, we shall just discuss the step-function 
potential case, leaving other more interesting potentials for the future. For this potential, all local 
densities are equal to their mean-field counterparts and, therefore, also the in-medium constants 
A = A(/r), P = P(/r) are expected to be spatially constant. Indeed, using the equilibrium (no 
medium) expressions Vo = Ub 0 (/to) = 7t 2 Bo(A/z/) and Eo = Eb q {h o) = 2tt 2 \uBo = po-Bo together 
with their in-medium values (see Eqs. (IV.4) and (IV.7)) 


Bb 0 (m) = e b 0 (h o)x (— + —) , Vb 0 (/t) = Vb 0 

1 VMo M J 
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we find the in-medium coupling constants 
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(V.19) 


(V.20) 


H* I 2 

Here we consider the liquid phase, i.e., /i > fio- In the gaseous phase, where the energy density 
as well as the particle density approach 0, the chemical potential has always the same value fio- 
Similarly, the in-medium skyrmion mass and its volume remain unchanged. 
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It should also be noticed that typical situations where such medium modified Skyrme La- 
grangians are considered are nucleons in atomic nuclei or in infinite nuclear matter at equilibrium. 
As we commented before, in this regime the BPS Skyrme model does not differ from its in-vacuum 
version. To make predictions for this regime, we have to include semiclassical corrections or extend 
the model to the near-BPS one. One the other hand, usually the coupling functions in the in¬ 
medium generalized Skyrme model are assumed to be spatially constant for infinite nuclear matter 
m ■ Here, the BPS Skyrme model can be of some help, at least for sufficiently flat potentials, 
because the original Skyrme model with its crystalline structure for large B definitely does not 
lead to a flat energy density. 

C. The BPS Skyrme model, a and oj mesons, and chiral symmetry 

Usually, in the Walecka model (or other, more general low energy effective models) the equilib¬ 
rium at nuclear saturation is the result of a precise balance between the repulsive forces induced 
by the ui meson and the attractive forces due to the a mesons (plus some small contributions of 
further mesons in more general models). It is part of the elegance of the BPS Skyrme model that, 
as a result of the BPS equation, it provides an exact cancelation between these forces, without any 
need for a fine-tuning of coupling constants, which leads to an exact equilibrium at the nuclear sat¬ 
uration density, and to exactly zero classical binding energies. Indeed, the static energy functional 
E = Eq + Eq = f d 3 x(se + £o); with eq = 7t 4 A 2 £>q and £q = v 2 U, consists of two terms which scale 
oppositely under Derrick scaling x —>• Ax. Concretely, Eq tends to expand the field configuration, 
inducing a repulsive force between different volume elements of the soliton, whereas the potential 
part Eq tends to collapse the configuration, corresponding to an attractive force. It is instructive 
to consider the resulting BPS equation for general pressure and for the axially symmetric ansatz 
(relevant, e.g., for neutron stars), which has a simple physical interpretation. Indeed, the 
equation reads 


(11.13 


£§{r) = £o(r) + P (V.21) 

or, in words, ’'the repulsive radial force per area equals the attractive radial force per area plus 
the pressure”. For P = 0, it just expresses the exact balance between repulsive and attractive 
forces at nuclear saturation, whereas in the limit of large pressure it shows that the repulsive force 
dominates, explaining the stiff character of the equation of state in that limit. 

As a result of the above, it is, therefore, reasonable to relate the two terms to the c o and a 
meson, respectively. Schematically we may write 

Cbps = - u 2 U = £ U (U) + C a (U). (V.22) 

Here, only the pionic degrees of freedom (chiral SU( 2) fields) are present explicitly, whereas the 
(or at least some) effects of the uj and a mesons are related to specific terms in the action. In other 
words, the u and cr mesons are hidden in the (nonlinear) Skyrme model action, and their effects are 
unravelled by studying the properties of particular solutions, similarly to the baryons themselves, 
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which, too, are absent in the action and become visible only on the level of (solitonic) solutions as 
coherent superpositions of pion fields. In other words, both baryons (and atomic nuclei) and the 
u and <7 mesons are realized in the model as emergent objects in a nonlinear pionic fluid. 

Let us remark that the exact balance between attractive and repulsive forces in the BPS Skyrme 
model will be destroyed if we add the perturbative part to the action, that is, extend the model to 
the near-BPS Skyrme model (or include the semiclassical corrections). However, one can control 
this transition using the small e parameter in the full action. 

As shown already, the sextic term leading to repulsion is responsible for the equivalence of the 
BPS Skyrme model and the Walecka model at high density (pressure), which allowed us to estab¬ 
lish a precise relation between the Skyrme model parameters and parameters of the u meson. This 
quantitative relation is possible because the sextic term and the omega meson give the leading 
contribution at large density in the two models. For the potential term, though, a quantitative 
relation to the a meson is not so obvious, because their contribution is subleading at large densi¬ 
ties, whereas at lower densities the a meson contributions in the Walecka model mix with other 
contributions (e.g., in-medium fermion contributions). 

Finally, let us recall that the potential term explicitly breaks chiral symmetry in the Skyrme 
model EFT, therefore it should also be related to the (spontaneous or dynamical) chiral symmetry 
breaking in the underlying fundamental theory, i.e., QCD, and should, therefore, depend on the 
corresponding order parameter (the quark condensate ( qq ); we remark that in the Walecka model 
there exists a linear relation between the nucleon condensate ( ijj'ip ) and the a meson VEV). The 
BPS Skyrme model without potential, on the other hand, which should correspond to the case 
without chiral symmetry breaking, behaves completely differently at low densities (e.g., there is 
no nuclear saturation), and only in the limit of infinite density (infinite pressure) the BPS Skyrme 
model approaches the behavior of the pure sextic model without potential. This implies that for 
finite density (finite baryon chemical potential) chiral symmetry remains broken, and no phase 
transition to a chirally symmetric phase occurs. 

Let us also notice that other types of phase transitions are absent, too, in the BPS Skyrme model. 
There is, e.g., no phase transition of a topological type - skyrmions always remain skyrmions, 
and no fractional topological state, for example half-skyrmion state, is created when the pressure 
is increased [36], [38], [39]. The creation of such phases probably requires the inclusion of the 
perturbative (i.e., the non-BPS) part of the full near-BPS action. The details of the transition to 
these new topological phases should, then, be related to the mutual strength of BPS and non-BPS 
parts of the full effective theory. It is interesting to note, however, that the introduction of an 
external magnetic field may lead to topological phase transitions, at least in 2+1 dimensions m- 
Also a phase transition to, e.g., quark matter cannot be described within the BPS Skyrme model 
alone. In other words, matter described by the BPS Skyrme model is always in a hadronic phase, 
and the only phase transition is the one between a gaseous hadronic phase below nuclear saturation 
and a liquid hadronic phase at and above nuclear saturation. It is interesting to note that this is 
precisely equivalent to the liquid-gas phase transition of nuclear matter, thus exactly reproducing 
the conjectured phase diagram of QCD at zero temperature for not too high values of the baryon 
chemical potential (close to nuclear saturation), see, e.g., [H] . 
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VI. BPS SKYRME MODEL AS A PERFECT FLUID 


A crucial property of the BPS Skyrme model was the fact that it has the energy momentum 
tensor of a perfect fluid, and that the static energy functional is invariant under SDiff transforma¬ 
tions on physical space. Here we want to show briefly that the relation goes much further and that, 
at least formally, the action of the BPS Skyrme model is equivalent to the action of a field theoretic 
description of perfect fluids in an Eulerian formulation [32]. There exist two main formulations of 
fluid mechanics, namely the Lagrangian formulation, where the dynamical variables are given by 
the particle trajectories (for finitely many particles) or by the fluid element trajectories (contin¬ 
uum limit), and the Eulerian formulation, where the dynamical variables have a more collective 
character and are provided by the (particle or mass) density p and by the fluid velocity v (see, 
e.g., pEHj). For finitely many particles, the d.o.f. in the Lagrangian formulation are the N particle 
trajectories X n (t), (n = 1,..., N), and the corresponding mass density in the Eulerian formulation 
is p(t,x) = mJ2n^ 3 \X n (t) — x). In the continuum limit, the discrete particle label n is replaced 
by three continuous labels y a , a = 1,2,3 required to label all fluid elements in three-dimensional 
space, and the corresponding dynamical variables are X(t,y) in the Lagrangian formulation and 

p(t,x) = po j d 3 y 3 (3) (x(t,y) - x) , 

v(t,x) = p~ 3 j where j = Po f d 3 y XS^ (^X(t,y) — of) (VI.1) 


in the Eulerian formulation. Here, the fluid element labels y a may be identified with the comoving 
coordinates of the fluid. For later convenience, we prefer to interpret p as a particle number density 
(not a mass density). In addition, we prefer to include the dimensions of the y a into po, such that 
the y a are dimensionless, which makes po dimensionless, too. We further remark that actions in the 
Lagrangian formulation, based on the X l (t,y a ) require integrations over y a , i.e., not over physical 
space. Actions for the Eulerian formulation, on the other hand, include integrals over physical space 
x l . The disadvantage of an action principle in the Eulerian formulation based on the dynamical 
fields p and v is that the constraints required by hydrodynamics (particle number conservation, 
etc.) require the introduction of Lagrange multipliers, which complicates the analysis. Recently, 
however, a field theoretic version of the Eulerian formulation of fluid dynamics gained support [32], 
where the constraints are satisfied identically, evading thereby the necessity of Lagrange multipliers. 
In this field theoretic version, the comoving coordinates y a are promoted to the dynamical variables 
of the theory. Indeed, for a regularly flowing fluid the flow function x l = X l [t,y a ) has an inverse 
y a = cj) a {t,x l ) which allows to express the density like 


p{t,x) 


poD where 


D = det 



(VI-2) 


Here we assumed an Euclidean target space and cartesian coordinates with a volume form equal 
to one, but it is convenient to allow for non-cartesian coordinates (and, eventually, for a curved 
target space) with volume form f l[(j> a )d(l) 1 d(j) 2 d(J) 3 , leading to 

d£\ 

dx l ) 


D = £l(4> a ) det 


(VI-3) 
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The determinant D has the further expression D 2 = Q 2 det(<9 i<f) a d l <p b ) which is useful, because it 
immediately allows for the relativistic generalization 

' d$ a dcj) b ' 


D 2 = Q 2 det 


dx p dx. 


(VI-4) 


Finally, there exists a third expression for D, 6 D 2 = AAVVV in terms of the particle number current 


A^ = ^P°e abc dA a dA h dA c - 


(VI-5) 


It is now easy to find the relativistic generalizations of v and p. The velocity is replaced by the 
four-velocity u A The cf) a are the comoving coordinates which do not change along the flow, which 
implies u^d^<j) a = 0. u p must, therefore, be proportional to the particle number current Af p , and 
the condition u p u p = 1 leads to 
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m _ 


Af p 
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\/(>D 


ne^e abc d u rd P ^ b d^ c , 


(VI-6) 


and the particle number density is defined by 


^ = pu" => p= pjM, = V6D = Ven(det(^^y. (VI.7) 

Observe that particle number conservation d^A = 0 (covariant conservation V^A= 0 in the 
general-relativistic case) is now an identity and does not require Lagrange multipliers. 

In this setting, a perfect fluid action is defined by choosing a Lagrange density F depending on 
(j) a , d p (j) a , etc. The simplest choice assumes that F depends on the target space variables only via the 
scalar p, but more general perfect fluids may depend on further thermodynamic variables Jia( 4 >a )- 
Here we shall permit at most one further thermodynamical variable h(cj> a ), i.e., F = F(p,h(cj) a )). 
The action just reads (we momentarily assume a general metric, for convenience) 


S = j d 4 xyAg\F(p,h) 


(VI-8) 


and leads to the energy-momentum tensor of a perfect fluid 


where 


T pa = -2\g \~3 -- S = (p + e)u p u a - pg pa 

bgpo 


_,. ,. de 

s = -F(p,h), p = p—-e. 


(VI.9) 


(VI. 10) 


A particularly simple case occurs for F = F(p). Then, both e = e(p) and p = p(p) are functions of 
the particle density p only, and, upon eliminating p , an equation of state e = e(p) may be found. 
Such fluids are called ”barotropic”. In the general case F = F(p , h), both e and p are functions of 
the two thermodynamic variables p and h. A natural choice is h = s(4> a ), where s is the entropy 
per particle. Observe that the entropy current S p = sAf p is conserved identically (d p S p = 0 or 
\7 P S P = 0, respectively). 









24 


It is now quite obvious how to relate the perfect fluid field theory sketched above to the BPS 


Skyrrne model. We just have to identify the Skyrme field U (the fields £, 0, <f>, see (11.12)) with the 


three scalar functions (f> a of the fluid. This identification is formal (only possible locally), because 
the cj> a take values in M 3 (or a subspace thereof), whereas the Skyrme field takes values in SU(2) 
(or, equivalently, in 8 3 ). If we further assume that the volume form 12 is - at least locally - the 
volume form on S 3 , 12 = sin 2 £ sin Qd^dQdQ, then the baryon density ps of the Skyrme model 
formally coincides with the particle number density p of the perfect fluid, and the baryon current 
Bv coincides with the particle number current AAh Globally, the two currents are, of course, 
different. The topology of § 3 guarantees, e.g., that the resulting baryon number B is always an 
integer, which is not true for the particle number current. If we accept this formal analogy, then 
the lagrangian of the BPS Skyrme model is related to the fluid lagrangian 


F = -A0tV - v l U{<\> a ) 


(VI-11) 


where p is the particle density and the potential U corresponds to a further ’’thermodynamical 
variable”. While the identification of p with the baryon density ps is obvious, it is not so clear 
what thermodynamical variable should correspond to the potential. The identification of U with 
the entropy per particle is not plausible, because the nuclear matter which the BPS Skyrme model 
is supposed to describe is essentially at zero temperature. Finally, the two cases of the BPS Skyrme 
model with the step function potential and without potential lead to lagrangians F(p) which only 
depend on p, and, therefore, correspond to barotropic fluids. BPS Skyrme models with genuine, 
(7-dependent potentials, on the other hand, correspond to non-barotropic (i.e., baroclinic) fluids. 


VII. CONCLUSIONS 

The main result of the present paper is the introduction and analytical description of the baryon 
chemical potential for the BPS Skyrme model, which is one necessary ingredient for a full under¬ 
standing of the thermodynamical properties of skyrmions as nuclear matter at zero temperature. 
We found the especially simple result that the baryon chemical potential is just the baryon charge 
density multiplied by a constant. As for other global quantities in the BPS Skyrme model (the total 
energy and volume), also the mean-field baryon chemical potential can be analytically obtained 
(as a function of the pressure P) by a certain integral (average) over the target space. In other 
words, the BPS Skyrme model realises to the very extreme the concept of a geometric model of 
matter. Indeed, all global quantities (energy, volume, mean-field energy density, mean-field chem¬ 
ical potential) can be computed without the knowledge of particular solutions. We also confirmed 
that the baryon chemical potential obeys all required thermodynamical relations. 

The complete analysis of the thermodynamics of the BPS Skyrme model at T = 0 also allowed 
to conclude that the model is equivalent to the Walecka model (as well as a vector interaction 
enhanced bag model)) in the high density regime. At densities close to the saturation point, 
however, the behavior of the BPS Skyrme model is more involved. Depending on the form of 
the potential or, more precisely, on how it approaches the vacuum, one may have a kind of bag 
model (with hadronic degrees of freedom) behavior with a non-zero energy density at zero pressure 
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or, instead, a zero energy density onset. The first case occurs for potentials whose near-vacuum 
dependence is U ~ a < 6 . This case is still qualitatively similar to the Walecka model in that 
both models show nuclear saturation at the saturation density. The main difference is that in the 
Walecka model there is a region of negative pressure (long-range attractive forces) below but close 
to nuclear saturation density, whereas pressure is exactly zero below nuclear saturation in the BPS 
Skyrme model, based on its classical soliton solutions. This behavior should, however, change once 
further terms of the near-BPS Skyrme model are included. It is, for instance, known that the 
standard non-linear sigma model term £2 = A 2 L a L a induces attractive long-range forces between 
skyrmions in some (attractive) spin-isospin channels [33]. Also the inclusion of binding energies 
due to quantum corrections (semi-classical quantization) of the skyrmion energies should change 
this behavior in a similar fashion. 

It is important to notice that the equivalence with the Walecka and related models is obtained 
in the mean-field limit. The BPS Skyrme model, on the other hand, also provides in a natural way 
a description beyond mean-field theory, that is, on the level of local, field theoretical quantities 
(energy density, baryon number density). In the asymptotic large density regime, where the exact 
equivalence is established, the local field theoretical computations agree with their mean-field 
approximations. In a lower density regime, however, where fi and /2 are different, we think that 
the full field-theoretic non-mean-field quantities (baryon chemical potential, energy density, etc.) 
should be used for a more precise description. 

Furthermore, we found some evidence that the BPS Skyrme model can be interpreted as an 
w-meson dominated model of nucleons, where the cu-mesons are hidden in the form of the action 
(or, more precisely, in the derivative term used in the action - baryon current squared) rather than 
included as an effective low energy field. We remind the reader that, due to the absence of the 
quadratic term £2 = A 2 there are no propagating pions and, therefore, no forces mediated by 
pions in the BPS Skyrme model. This simply means that, as stated repeatedly, the BPS Skyrme 
model by itself cannot be considered a complete low-energy effective theory for nuclear physics or 
strong interactions. A good candidate for such a low-energy effective theory is, in our opinion, the 
near-BPS Skyrme model 

£ = £bps + e (-£2 + £4 + A>) (VII. 1) 

(here £4 = A 4 [£ M , L u ] 2 is the Skyrme term, and £0 is a further potential), where e is assumed to be 
small in the sense that £bps provides the main contributions to soliton masses and is dominant in 
regions of sufficienty large baryon density, such that the unique properties of the BPS submodel are 
essentially preserved in this regime. Close to the vacuum, on the other hand, £2 always dominates 
over £ 6 , and long-range forces mediated by pions are, therefore, present in the near-BPS extension. 
It is of some interest to note that the BPS Skyrme model does reproduce the forces related to the 
u and a mesons - whose fields do not appear in the action - while it does not include the pionic 
forces, although its action is expressed entirely in terms of pion fields. 

As already mentioned, there are not many results concerning the baryon chemical potential 
in the Skyrme theory. One comparison, however, may be made with results found by means of 
the AdS/CFT correspondence and the Sakai-Sugimoto model |45| . which is a holographic (large 


26 


N c relevant) version of a Skyrme type model. It has been found that for the two flavour case 
asymptotically the baryon density behaves as for a free fermion gas ps oc /i 3 [36], which strongly 
differs form our result. In spite of that, there is a qualitative similarity between this holographic 
computation and our approach. In both models, the baryon density is always a function of the space 
coordinates at any finite value of the chemical potential. Such an inhomogeneous configuration 
tends to a homogeneous one only in the infinite density (chemical potential) limit [36]. On the 
other hand, it is possible to get an asymptotically linear relation between the chemical potential 
and baryon density in a holographic set-up. For this, one needs the usual Maxwell action instead 
of the Born-Infeld one m- Adding more scalar fields may, however, change this linear relation. 

There are several obvious directions in which the current investigations should be continued. 
First of all, the knowledge of the baryon chemical potential for the BPS Skyrme model is essential 
if one wants to apply the model for a complete description of neutron stars. In principle, the core 
of neutron stars, described by the BPS Skyrme model, may be surrounded by a skin with more 
usual matter with a known equation of state. The obvious condition for a transition from the 
dense hadronic phase (BPS Skyrme action) to the skin phase is the equivalence of the chemical 
potentials. This should lead to a modification of the mass-radius relation for the low massive stars 
with perhaps an appearance of a minimal neutron star mass [38]. This, however, will be modified 
already by the inclusion of the non-BPS part of the full near BPS Skyrme action. Indeed, as 
one approaches the outer region of a neutron star in the BPS Skyrme model, the matter Skyrme 
field tends to its vacuum value. But close to the vacuum, the perturbative terms in the chiral 
Lagrangian dominate over the BPS part. Hence, in this regime they cannot be omitted. This 
points towards another important issue which should be understood, namely the generalization of 
the thermodynamical description of the BPS Skyrme theory to its near-BPS extension, as well as 
its application to nuclear phenomenology. Let us notice that, except for the step-function potential 
case, one should use the baryon chemical potential rather than its mean-field approximation. In 
fact, for non-flat potentials it is known that local quantities, especially for heavy neutron stars, 
change a lot if computed in the mean-field limit (compare pressures and energy densities inside 
neutron stars in [29]). 

Another important quantity which has to be understood is the isospin chemical potential (there 
already exist some proposals on how to treat this issue in the Skyrme framework [49] ). This would 
allow for a complete description of skyrmionic nuclear matter at zero temperature. 
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